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We report experimental observations of an undulational instability of myelin figures. Motivated 
by this, we examine theoretically the deformation and possible instability of concentric, cylindrical, 
multi-lamellar membrane structures. Under conditions of osmotic stress (swelling or dehydration), 

we find a stable, deformed state in which the layer deformation is given by 5R (x rV-^^/(''-^) , where 
Ba is the area compression modulus, B is the inter-layer compression modulus, and h is the repeat 
distance of layers. Also, above a finite threshold of dehydration (or osmotic stress), we find that 
the system becomes unstable to undulations, first with a characteristic wavelength of order y^^do, 
where ^ is the standard smectic penetration depth and do is the thickness of dehydrated region. 

PACS numbers: PACS numbers: 87.22. Bt, 64.60.-1, 82.70.-y, 68.15.-fe 



The swelling of bulk phospholipid by water at a tem- 
perature higher than the chain melting temperature can 
lead to the growth of multilayer hpid tubules from the 
interface of bulk lipid and water j^, ||, ^,0, 1). These 
tubules are of mesoscopic dimensions with a diameter 
about 20-40 fim. They have the symmetry of a smectic- 
A liquid crystal and a strong optical anisotropy. Both 
electron microscopy Q and X-ray diffraction show a con- 
centric multilayer structure of these lipid tubes with a 
layer spacing about 60A |^]. Due to the similarity of 
their structural features to those of nerve myelin sheaths, 
these tubules are often referred as "myelin tubes" or 
"myelin figures". The morphological features of myelin 
tubes have been classified into two steps according to the 
time of growth During the first step, simple tubes 
grow into bulk water with an initial growth rate in length 
about 1 fim/sec. The diameter and number of layers of 
tubes appear to remain constant during this process and 
the growth rate is inversely proportional to the square 
root of growth time ^. The growth almost stops in 
the second step which is characterized by the formation 
of complicated morphologies of myelin tubes, such as he- 
lical and coiling forms, apparently in order to maximize 
their inter-membrane attraction g. 

In our experiments, we have observed an instability 
of myelin tubes under dehydration by isolating a single 
myelin tube from others (during the first step, as shown 
in Fig. 1) and allowing the bulk water to evaporate in a 
controlled open chamber. During this dehydration step, 
periodic bumps with a wavelength about 1 /im are ob- 
served on the surface of myelin tubes as shown in Fig. 
2(a). As water further evaporates, these bumps grow 
into arms, as shown in Figs. 2(b) and 2(c), and a simi- 
lar instability can occur for these long arms. As far as 
we know, this is the first reported observation of the in- 
stability of myelin tubes under dehydration. After the 
bulk water dries out, the tubular structure of myelins 



disappear and a mosaic structure is observed. Although 
the full structural change of myelin tubes under dehydra- 
tion is very complicated and highly nonlinear, the initial 
instability of myelin tubes under dehydration poses an 
interesting deswelling of layered structures in cylindrical 
geometry, which we examine theoretically. 

Here, we suggest an explanation for the initial morpho- 
logical change of myelin tubes observed in the dehydra- 
tion step based on energetic arguments. We first demon- 
strate and examine a stable configuration of myelin un- 
der conditions of osmotic stress, in which layer bending 
does not occur. We then examine the limit of stability of 
this to undulations, which relieve in-plane compression. 
We predict the appearance of undulations with a finite 
wavelength that can be of order micrometers. We then 
speculate on the kinetic aspects of the deswelling pro- 
cess. We note that, unlike the pearling instability pO[ , 
the compression and bending energies already favor finite 
wavelength undulations in this case, and the appearance 
of such an instability need not be due to hydrodynamic 
effects. These, we argue will not significantly change our 
results for slow drying rates. 

For a uni-lamellar tube, dehydration tends to decrease 
its radius but increase its length. Surface undulation 
is disfavored because it increases tube's bending energy. 
For a multi-lamellar tube, in addition to the bending en- 
ergy of each layer, inter-layer and in-plane compression 
energies need to be considered. It can be shown that, in 
the limit of long cylinders, a change in length becomes en- 
ergetically prohibitive. We show below that undulation 
of dehydrated layers occurs only for finite dehydration of 
a multi-lamellar tube in the presence of an undisturbed 
core, while in-plane compression occurs (without undula- 
tions) below this threshold. Although a one-dimensional 
undulation along the long axis of myelin tubes is reported 
in this paper, we cannot exclude the possibility of a gen- 
eral two-dimensional undulation along the surface of the 
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tube. 

We begin our discussion by considering the dehydra- 
tion of a myehn tube in bulk water. The length L of 
the myelin tube is assumed to be much larger than its 
radius i?o- As water evaporates, the osmotic pressure 
changes due to the increased impurity concentration in 
water. In response to the drop in chemical potential out- 
side, water molecules in the outer layers of the myelin 
tube will diffuse out to the bulk and the myelin tube 
will deform. This process is relatively fast, as the mem- 
branes are permeable to water. Migration of lipids will 
be ignored here, since this is a much slower process due 
to their low solubility in water. The removal of water 
thus results in a radial deformation 5R < of a cylin- 
der of initial radius R. Through second order in this 
deformation, the area change of this cylinder is given by 
M = 27r / dz {6R + R\V6R\^/2). Here, the gradient op- 
erator is taken to act in the plane of the membrane. 

We shall find below that there is an initial uniform (in 
z) deformation SR that, beyond a finite threshold, drives 
an undulatory instability. We take this undulatory con- 
tribution SR' to be sinusoidal. Thus, in order to deter- 
mine the onset of this instability, we shall keep terms 
only through second order in SR' . With this in mind, we 
take the area strain to be e = SR/R + \\/SR'\'^/2. The 
first term in e is the area deformation due to the average 
radius change in a cylinder, while the second term is due 
to layer undulation This strain is actually averaged 
over one period along the cylinder axis. For fluid mem- 
branes, it would not be correct to integrate the square 
of the local area strain, as molecules are free to migrate 
within each layer. Rather, the integrated square of e 
above properly accounts for the energy associated with 
area compression, through second order in SR' . 

We parametrize the structure by the radial position r 
before dehydration, and the position z along the axis Jill. 
Near the onset of the instability, the displacement field 
SR{r, z) = SR + SR' is expected to consist of a displace- 
ment SR{r) that is uniform along the cylinder axis plus 
a sinusoidally varying displacement SR'{r,z) oc cos(gz). 
For such a deformation, the associated energy density in 
our model can be expressed as 



+ iK{V'SR')\ 



(1) 



where Ba is the in-plane compression modulus, B is the 
inter-layer compression modulus, K is the bulk bending 
rigidity, and h is the repeat distance of layers before de- 
hydration. 

Before discussing possible undulational instabilities of 
the myelin tubes, we first consider deswelling of the struc- 
ture in the absence of undulations {i.e., SR' — 0). Here, 
we ignore the bending contributions. (Note the small cur- 
vatures, with radii of order several micrometers, involved 
in especially the outer layers of the myelin, where the 
deswelling is greatest.) In this case, the shift in chemical 



potential of water molecules in a layer at radius r' due 
to layer deformation can be calculated as 
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where F = J (2TTrf)dr is the free energy per unit length, 
N(r') represents the number of water molecules per unit 
length in a layer at radius r', and vo is the molec- 
ular volume. In Eq. (^), we have used the relation 
27rr [(5i?(r)] = VQSN{r')6{r — r') resulting from simple 
geometry, where the Heaviside step function 9{r — r') en- 
sures that 5R{r) can only be changed by changing N at 
r' < r. Uniformity of the chemical potential in equilib- 
rium (at least over time scales for which lipid transport 
is negligible) requires the integrand in Eq. to vanish. 
The solution to this differential equation is 



SRo^r", 



(3) 



where a = y/ Ba/ (Bh). This result applies for both 
swelling {SR > 0) and dehydration {SR < 0), where the 
amplitude of the displacement field is determined by the 
value of the chemical potential outside the cylinder. The 
result is valid near the surface where bending contribu- 
tions are negligible, and is sensible in that for large Ba 
there will be little penetration of the dehydration profile 
in from the surface, while for large B, the profile is more 
uniform toward the outside of the cylinder. 

Upon changing the osmotic pressure in the bulk water 
a dehydration profile will extend throughout the cylin- 
drical structure. The initial stage is a smooth deswelling 
with only a change in the radius of each layer as calcu- 
lated above. This also involves bending energy, as the 
dehydration profile extends deeper and becomes more 
pronounced. For a finite degree of deswelling, however, 
the layers can bend to relieve the intra- and inter- layer 
compression energies. Hence, we now show that incor- 
porating the bending degrees of freedom leads to an 
instability to a surface undulation. For simplicity, we 
assume a single mode approximation and a linear de- 
formation profile that extends to a depth do below the 
outer surface Rq of the myelin, i.e., 5R{r = Rq — x) — 
(1 — x/do) (a -I- a'cos(g2:)), where a and a' are the uni- 
form and undulatory deformation amplitudes, x is the 
depth below the outermost layer, and q is the wave num- 
ber (as shown in Fig. 3). This approximates the state 
of the system at an instant in time during dehydration, 
rather than an equilibrium profile. 

In this case, the modulated part of the total free energy 
per unit length can then be calculated from Eq. Imj as 
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Minimizing Eq. (|[) with respect to the wave num- 
ber (7, the optimal wave number is given by (9*)^ — 
3BAa/8hKRo. Such an undulation is favored only when 
Ai^,„nH < 0. This results in a threshold 
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below which the system is stable. The corresponding 
reduction in volume of the cylinder is then given by 
AVc — 27rRQLac- At the onset of instability, we find 
that the wavelength of the instability is 

A ~ (7) 

which depends only on material parameters and on the 
depth of the deswelling profile. Note that any finite 
range dehydration profile would lead to a similar result, 
with different prefactors. For lipid bilayers, the stan- 
dard smectic penetration depth ^ = -y/ K/B is about 0.01 
/im. Therefore the undulation wavelength is of order 1 
/im if do is a few microns. We have considered a one- 
dimensional undulation along the long axis of the myelin 
tube. Additional undulation along the ip direction might 
exist. However, our calculation shows that such an un- 
dulation is disfavored if the undulation amplitude and 
wavelength are isotropic in both directions. We also note 
that a change in layer spacing due to surface area com- 
pression is much smaller than that due to dehydration 
(their ratio is less than do/Ro) and can be neglected. 

The discussion thus far has been based on energetic 
considerations. A dynamic calculation would incorporate 
the energetic driving force for the instability and dissipa- 
tive damping, and conservation of water in a given layer 
during the instability should shift the wavelength A to a 
smaller value jl^. Moreover, there is the possibility that 
large effluence of water induced by a severe dehydration 
could induce a dynamic, rather than quasi-equilibrium, 
instability through non-linear convective couplings fl^ . 

In the early stage of dehydration the dominant pro- 
cess limiting diffusion is the compression energy between 
layers, and the associated diffusion equation for water 
molecules to cross bilayers near the surface r ~ Rq is 
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where (/> is a coarse-grained composition variable repre- 
senting the relative change in water composition (ratio 
of interlayer spacing after dehydration to spacing be- 
fore dehyeration) , and Z? is a mobility proportional to 
the membrane permeability. As dehydration proceeds 
the intra-layer compression elasticity Ba will eventually 
slow this process. In our calculation above we have ap- 
proximated the profile </)(r, t) by a piecewise linear profile 
of depth do, but other bounded profiles lead to similar 
results. At a time t after the external osmotic pres- 
sure has been quenched, the penetration of the dehy- 
dration layer will thus scale as d^it) ~ ^BDt. Using 



the piecewise-linear approximation above for 5R, the to- 
tal volume change of the outermost cylinder is given by 
AV = nLRodo (1 — 0o)- The change in layer spacing of 
the outermost layer is dR{r) — 5R{r — h) ~ h{(j)o — 1). 

Our instability condition, Eq. (||), assumes a given 
profile and dehydration depth. Assuming further the 
local equilibrium assumption, that the chemical poten- 
tial of the outermost layer is equal to the chemical po- 
tential ^ext specified by the bulk, leads to the relation 
fiext — —voB(l — (po). (Corrections to this involve terms 
that go as l/i?o and l/i?o from Eq. (||)). A typical sce- 
nario would involve sealing the chamber after a rapid 
dehydration to maintain, for a large water-to-lipid ra- 
tio, a constant osmotic pressure 11 — ^ext/vo outside, 
and hence the relation {1 — cpo) ^ H/B at the outermost 
layer. In this case the critical size of the dehydrated re- 
gion is do — yJ\Qh£,RQB^I{BAa). Equivalents, the crit- 
ical time scale for the instability to occur after the dehy- 
dration is r ~ h£^B Rq / {D B A^) , which increases linearly 
with Rq. The initial undulation wavelength is propor- 
tional to -Rq^^. 

One could envision another extreme of imposing a con- 
stant dehydration flux j by, for example, exposing the 
chamber to an unsaturated vapor that continually in- 
creases the impurity concentration Cimp in the bulk wa- 
ter, hence increasing the osmotic pressure and the driv- 
ing force for dehydration. In this case the details of 
the timescale necessary for the instability would depend 
on comparing the rate of change of (f>o to the time for 
diffusion out of the multilamellar state, although the 
relation between outer composition and overall volume 
change, AVc = 2TrRoLac, would still hold. For the out- 
ermost layer, the stress of layers must balance the os- 
motic pressure difference to the bulk water. In this case 
(fiQ (t) ^ I — jt, for a constant flux j ^ kTcimp/ B where 
kT is the thermal energy, and the critical value of (pQ 
is ~ 1 - 10BhS,Ro/ {BAdl) . The delay time before 

undulation is r ~ \/h^Ro/ {BaDJ) oc ■\/i?o/i- The ini- 

1/8 

tial undulation wavelength is proportional to {Rq/ j) 
At a higher drying rate, the thickness of the dehydrated 
region is smaller at the onset of instability which corre- 
sponds to a smaller onset undulation wavelength. There- 
fore, in an experiment with controllable humidity, the 
undulation wavelength can be varied by controlling the 
relative humidity. Nevertheless, we do not have quanti- 
tative measurements yet. 

To conclude, we have studied the deformation of myelin 
tubes under dehydration by curvature elasticity. The 
cylindrical geometry is stable for infinitesimal dehydra- 
tion, and the surface area of dehydrated layers is reduced 
due to an in-plane compression. If an equilibrium is 
achieved at this stage, the equilibrium radius deformation 
profile can be described by 5R (x ^^/(hB) ^ p^^, g^ite 
dehydration, the cylindrical geometry becomes unstable 
and undulation occurs to reduce the in-plane compression 
energy. The critical volume change for this instability is 
estimated to be AVc ^ SOh^LRlB/{doBA)- The undu- 
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FIG. 1: A low resolution view of myelin tubes during the 
first step. The diameter of tubes is about 20-40 microns. An 
isolated myelin tube can be obtained by introducing a water 
flow across the chamber to disintegrate the bulk lipid. 



lation wavelength is predicted to be about 5\/^3o and is 
comparable to the observed undulation wavelength. 
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FIG. 2: A series of photographs of a myehn tube showing the 
instabihty under dehydration at different times. The wave- 
length of undulationis about 1 fim. At later times, the am- 
plitude of undulation increases and a bump becomes an arm. 





FIG. 3: Schematic representation of the undulation profile of 
a myelin tube: a cross-sectional view showing different layers. 



